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Sialidase catalyzes the removal of a terminal sialic acid from 
glycoconjugates and plays a pivotal role in nutrition, cellular 
interactions and pathogenesis mediating various infectious 
diseases including cholera, influenza and sepsis. An array of 
antiviral sialidase agents have been developed and are 
commercially available, such as zanamivir and oseltamivir 
for treating influenza. However, the development of bacterial 
sialidase inhibitors has been much less successful. Here, 
natural polyphenolic geranylated flavonoids which show 
significant inhibitory effects against Cp-NanI, a sialidase from 
Clostridium perfringens, are reported. This bacterium causes 
various gastrointestinal diseases. The crystal structure of the 
Cp-NanI catalytic domain in complex with the best inhibitor, 
diplacone, is also presented. This structure explains how 
diplacone generates a stable enzyme-inhibitor complex. These 
results provide a structural framework for understanding the 
interaction between sialidase and natural flavonoids, which are 
promising scaffolds on which to discover new anti-sialidase 
agents. 
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1. Introduction 

Viral and bacterial sialidases have been focused on as drug 
targets for the treatment of human infections (Soong et al, 
2006; Memoh et al, 2008). Viral sialidases are critically 
required for viral propagation because they are required for 
release from the host cell (Kim et al., 2013). Drug discovery of 
inhibitors targeting sialidase has focused on either synthetic 
compounds through random screening or substrate (A^-acetyl- 
neuraminic acid) mimics. However, these studies have mainly 
been targeted against siahdase from the influenza virus. There 
are several inhibitors that mimic the substrate, including 
zanamivir (Relenza), oseltamivir (Tamiflu; a cyclohexene 
derivative), peramivir (a cyclopentane-based inhibitor) and 
A-315675 (a pyrrolidine-based inhibitor; Kim et al, 2013). All 
of these compounds are based on 2-deoxy-2,3-dehydro-A^- 
acetylneuraminic acid (Neu5Ac2en), a putative transition- 
state analogue. At present, the existence of drug-resistant 
strains against zanamivir (Gubareva et al., 1998), oseltamivir 
(Gubareva et al, 2001) and other antivirals (Tambic Andra- 
sevic, 2004) have spurred on the need to discover innovative 
antiviral agents. 

On the other hand, pathogenic bacteria having resistance 
to pre-existing drugs occasionally cause fatal infectious 
diseases. The need to combat these superbugs mandates 
the development of new drugs overcoming antibacterial 
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resistance. Various pathogens such as Clostridium perfringens, 
Pseudomonas aeruginosa and Streptococcus pneumoniae 
express siahdases. These enzymes specifically catalyze the 
hydrolysis of terminal sialic acids from glycoconjugates 
(Shinya et al , 2006) and play pivotal roles in the pathogenesis 
of a number of microbial diseases, including cholera, enter- 
otoxaemia, gas gangrene and pneumonia (Corfield, 1992; 
Rood, 1998). Therefore, bacterial siahdases have recently 
emerged as a prominent target for the treatment of bacterial 
infections (Nguyen et al, 2010; Woo et al, 2011). The latest 
studies on sequelae of sepsis and septic shock have focused on 
understanding the roles of bacterial pathogenic factors 
(Rittirsch et al, 2008). Moreover, studies using a sialidase 
deletion mutant from P. aeruginosa showed that bacterial 
siaHdase acts in the initial stages of pulmonary infection by 
targeting glycoconjugates and biofilm production (Soong etal, 
2006). 

Since siahdase inhibitors are substrate analogues and typi- 
cally exhibit competitive kinetics, it is interesting to explore 
natural non-substrate mimics such as flavonoids, (oligo)- 
stilbenes, coumarins and diarylheptanoids (Grienke et al, 
2012). Flavonoids, which are polyphenolic compounds, are 
widespread in the plant kingdom and are known to show 
antibacterial and antiviral effects (Liu et al, 2008; Jeong et al, 
2009; Grienke et al, 2012). Natural flavonoids have also been 
reported to have latent anti-influenza and antibacterial effects 
against SARS (severe acute respiratory syndrome) virus, 
MRSA (methicillin-resistant Staphylococcus aureus) and 
Chagas disease mediated by a trypanosome (Arioka et al, 
2010; Cho et al, 2013; Navratilova et al, 2013). 

Recently, we identified the natural geranylated flavonoids 
diplacone, mimulone, 3^-O-methyldiplacone and A'-O- 
methyldiplacone from the fruit of Paulownia tomentosa 
(empress tree) and showed that they are potent inhibitors of 
both choHnesterase and butyrylcholinesterase (Cho et al, 
2012). Here, we report that geranylated flavonoids are novel 
inhibitors targeting bacterial sialidase. These compounds show 
striking competitive inhibitory effects against sialidase from 
C. perfringens (Cp-NanI). Importantly, C. perfringens is a 
human pathogen causing various gastrointestinal diseases such 
as enterotoxaemia, gas gangrene and peritonitis. 

In spite of reports of a number of structures of siahdases 
from bacteria (Crennell et al, 1993, 1994; Gaskell et al, 1995), 
viruses (Varghese et al, 1983; Burmeister et al, 1993; Crennell 
et al, 2000) and humans (Chavas et al, 2005), elucidation of 
the binding mode of bacterial siahdase inhibitors is limited 
owing to the absence of structural information (Taylor, 1996). 
In order to understand how natural geranylated flavonoids 
interact with siahdase, we determined the crystal structure of 
the Cp-NanI catalytic domain (Nanlco; residues 243-694) in 
complex with diplacone at 1.9 A resolution. The structural 
information on Cp-Nanlco and comparison with the models of 
human sialidases provide an insight into how the natural 
geranylated flavonoids can behave as antibiotics and led to the 
elucidation of a strategy for the design of advanced inhibitors 
for the treatment of infectious diseases caused by C. perfrin- 
gens. 



2. Methods 

2.1. Cloning^ expression and purification 

The DNA template for C. perfringens sialidase (Cp-NanI) 
has NCBI reference sequence WP_01 159033 1.1. For bacterial 
expression, the NanI catalytic domain without the lectin 
domain (NanlcD; residues 243-694) was amplified from 
C. perfringens genomic DNA by the polymerase chain reac- 
tion and cloned into the Ndel and Xhol restriction sites of the 
pET-23d expression vector (Novagen), which contains a direct 
C-terminal hexahistidine (Hisg) tag. Cp-NanlcD protein was 
expressed in Escherichia coli BL21 CodonPlus (DE3) cells. 
The transformed cells were cultured in Luria-Bertani (LB) 
medium containing 100 jig ml~^ ampicillin at 310 K. The 
bacterial cells were induced at 293 K for 20 h with 0.5 mM 
isopropyl y^-D-l-thiogalactopyranoside (Pharmacia) at an 
ODgoo of 0-6 and then harvested by centrifugation at 4000g for 
20 min. The harvested cells were resuspended in lysis buffer 
(50 mM sodium phosphate pH 7.0, 300 mM NaCl, 5 mM 
imidazole, 1 mM PMSF, 2 mM yS-mercaptoethanol) and lysed 
by sonication, after which the lysate was centrifuged at 14 OOOg 
for 1 h. The resultant supernatant was applied to an immobi- 
lized metal-affinity chromatography on nickel-nitrilotriacetic 
acid resin (Peptron) pre-equilibrated with lysis buffer (50 mM 
sodium phosphate pH 7.0, 300 mM NaCl, 20 mM imidazole). 
The column was then washed with ten bed volumes of wash 
buffer. The His6-tag fused protein bound to the column was 
eluted with elution buffer (50 mM sodium phosphate pH 7.0, 
300 mM NaCl, 300 mM imidazole). The samples were then 
purified by size-exclusion chromatography using a HiLoad 
16/60 Superdex 200 column (GE Healthcare Life Science) pre- 
equilibrated with gel-filtration buffer (25 mM CHES-HCl 
pH 9.5, 200 mM NaCl), after which the fractions containing 
Cp-NanlcD protein were collected. The protein was concen- 
trated to 21.8 mg ml~^ using an Amicon Ultra-15 30K (Milli- 
pore) and was stored at 193 K. 

2.2. Crystallization^ data collection and structure 
determination 

Crystallization trials were performed using the hanging- 
drop vapour-diffusion method. Crystals grew at 293 K in 2 ill 
drops comprised of equal volumes of protein solution and 
reservoir solution consisting of 20%(w/v) polyethylene glycol 
(PEG) 3350, 0.2 M ammonium sulfate. The crystals were 
soaked in soaking solution [0.1 M sodium cacodylate pH 6.5, 
20%(w/v) PEG 3350, 10%(>v/v) glycerol, 5 mM diplacone, 
5%(v/v) dimethyl sulfoxide] for 2 h at 293 K. The crystals were 
cryoprotected by transfer into cryoprotectant [20%(w/v) 
polyethylene glycol 4000, 10%(v/v) glycerol, 0.2 M ammonium 
sulfate] and flash-cooled in liquid nitrogen for data collection. 
A complete data set for Cp-NanlcD with diplacone at 1.9 A 
resolution was collected at 100 K using an ADSC Q315r 
detector on beamline 5C at Pohang Accelerator Laboratory 
(PAL) using an X-ray beam with a single wavelength 
(0.9795 A). The crystal belonged to space group P2i2i2i, with 
unit-cell parameters a = 69.1, b = 72.6, c = 97.1 A. Diffraction 
data were processed and scaled using the HKL-2000 suite 
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Table 1 

Data-collection and refinement statistics. 

Values in parentheses are for the highest resolution shell. 



Data collection 
X-ray source 
Wavelength (A) 
Resolution range (A) 
Space group 

Unit-cell parameters (A) 
No. of observed reflections 
No. of unique reflections 
Completeness (%) 



et (%) 



Mean I/cr(I) 
Multiplicity 
Refinement statistics 
Resolution range (A) 

^work/afreet (%) 

No. of atoms 

Protein 

Diplacone 

Ca^^ 

Water 
Average B factors (A^) 

Protein 

Diplacone 

Ca^^ 

Water 

R.m.s. deviations from ideal geometry 

Bond lengths (A) 

Bond angles (°) 
Ramachandran plot 

Most favoured regions (%) 

Allowed regions (%) 
PDB code 



Beamline 5C, PAL 
0.9795 

50-1.9 (1.93-1.90) 

a = 69.1, b = 72.6, c = 97.1 

195578 

39033 

99.1 (98.1) 

5.5 (9.6) 

32.7 (26.9) 

5.0 (5.1) 

40-1.9 
15.8/19.4 

3526 

31 

2 

574 

31.2 
50.1 
36.4 
41.7 

0.012 
1.30 

96.2 

3.8 

412e 



t Emerge = Em/ E. Ui(hkl) - {I{hkm/Y.hki Hihihkl), whcrc I{hkl) is the intensity of 
reflection hkl, Ylhu the sum over all reflections and E; is the sum over i measurements 
of reflection hkl. X i^^ork = Em/ ll^obsl - l^caicll/E^/ l^obsH; ^free is the R value 
calculated for 5% of the data set that was not included in the refinement. 



2.3. Preparation of geranylated flavonoids 

The extraction and purification of the natural geranylated 
flavonoids diplacone, mimulone, 3^-O-methyldiplacone and 
4^-O-methyldiplacone from P. tomentosa (empress tree) were 
performed as described previously (Cho et ai, 2012). 

2.4. Enzymatic inhibitory assay 

The sialidase inhibitory assay was performed as reported 
previously (Kim et aL, 2012; Cho et aL, 2013). Briefly, 15 jil 
GH33 sialidase solution (0.1 U ml~^) was pre-mixed with 15 |il 
sample solution at different concentrations in 510 |il 50 mM 
sodium acetate buffer pH 5.0 in a cuvette. 60 |il 0.125 mM 
4-methylumbelliferyl-Qf-D-A^-acetylneuraminic acid sodium 
salt hydrate (Sigma, catalogue No. M8639) in buffer at pH 5.0 
was then added to the mixture as a substrate to start the 
reaction at 310 K. 4-Methylumbelliferone was immediately 
quantified by fluorometric determination with a SpectraMax 
M^^ Multimode Reader (Molecular Devices, USA). The 
excitation wavelength was 365 nm and the emission wave- 
length was 450 nm. Enzyme activity was recorded over a range 
of pre-incubation times (Fig. la) or concentrations (Fig. lb). 
The data were analyzed using a nonlinear regression program 
(SigmaPlot; SPCC Inc., Chicago, Illinois, USA). The experi- 
mental enzyme activity was determined using the logistic 
curve represented by (1), based on a time-driven protocol with 
initial velocity. 



activity (%) = 100 x 



1 



[l + ([I]/IC5o)] 



(1) 



(Otwinowski & Minor, 1997). The structure of Cp-NanlcD 
with diplacone was determined using the molecular-replace- 
ment method with MOLREP (Vagin & Teplyakov, 2010) 
in the CCPA suite (Winn et al, 2011) using the apo crystal 
structure of Q^-NanlcD (PDB entry 2vk5; Newstead et al, 
2008) as the phasing model. Automated model building and 
refinement of NanlcD complexed with diplacone were 
performed with REFMAC5 (Murshudov et al, 2011) in the 
CCPA suite and phenix.refine (Adams et ai, 2010). phenix. 
ligandfit (Terwilliger et al, 2006) was used for automated 
ligand building of the diplacone in Q7-NanIcD, and manual 
model building was performed with Coot (Emsley & Cowtan, 
2004). Weak electron-density regions (residues 692-694) were 
excluded from the final structure. The Ramachandran statis- 
tics were calculated using MolProbity (Chen et al, 2010). 
96.2% of the residues were in the preferred regions and 3.8% 
of the residues were in allowed regions. A 2Fo — composite 
OMIT map of diplacone contoured at 1.0a was generated with 
FFT (Read & Schierbeek, 1988) in the CCPA suite. The atomic 
coordinates and structure factors of Q7-NanIcD with dipla- 
cone have been deposited in the Protein Data Bank (PDB) 
with accession code 412e. Data-coUection and refinement 
statistics are summarized in Table 1. 



2.5. Homology modelling 

Model structures of human sialidases (Neul, Neu3 and 
Neu4) were built using the MODELLER 9v7 software (Sah & 
Blundell, 1993). Human Neul-4 protein sequences were 
retrieved from NCBI database (Neul, NM_000434.2; Neu2, 
NM_005383.2; Neu3, NM_006656.5; Neu4, NM_001167599.1). 
On the basis of the Neu2 sequence and structure (PDB entry 
Ivcu; Chavas et al, 2005), multiple sequence alignment and 
structure-based alignment were performed using ClustalXl 
(Larkin et al, 2007) and MODELLER 9v7 alignld (Sah & 
Blundell, 1993; Supplementary Fig. Sl^). Regions from Neul, 
Neu3 and Neu4 that are not conserved in the Neu2 catalytic 
domain (Neul residues 1-64, Neu3 residues 1-11, 287-300 and 
315-326 and Neu4 residues 1-10, 284-336 and 355-373) were 
excluded. Modelling was carried out for the human Neul, 
Neu3 and Neu4 catalytic domains against the chosen Neu2 
template (PDB entry Ivcu) using MODELLER 9v7 model- 
single (Sah & Blundefl, 1993; Supplementary Fig. S2). Each of 
ten outputs of the modelled structures was ranked using an 
internal scoring function in MODELLER 9v7. The reliability 
of the homology modelling was assessed by calculating the 
root-mean-square deviation (r.m.s.d.) using PyMOL v.l. 5.0.4 

^ Supporting information has been deposited in the lUCr electronic archive 
(Reference: MV5098). 
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(DeLano, 2004; Supplementary Table SI), evaluation of the 
Z-score using ProSA (Wiederstein & Sippl, 2007) and analysis 
of the q) and ^jr torsion angles using PRO CHECK (Laskowski 
et al, 1993; Supplementary Table S2). 

3. Results 

3.1. Kinetic properties of the natural geranylated Cp-NanI 
inhibitors 

In order to demonstrate that flavonoid-based natural 
compounds extracted from P. tomentosa act as inhibitors of 
Cp-NanI, a comparative activity screen was performed. This 
method relies on substrate hydrolysis, leading to an increase in 
fluorescence that can be read out in real time. Interestingly, a 
geranyl (3,7-dimethylocta-2,6-dienyl) group at the C6 position 
within the flavonoid was found to be crucial for both inhibi- 



tory potency and kinetic mechanism. For example, diplacone 
showed a 150-fold higher efficacy (IC50 = 0.11 \iM) than 
eriodictyol, a parent flavonoid lacking a geranyl group (IC50 = 
16.4 |iM) (Table 2). As shown in Fig. 1, diplacone showed 
competitive kinetics, with a value of 0.04 |iM. Diplacone 
also acted as a time-dependent inhibitor. 

Slow-binding inhibition mechanisms can be investigated by 
pre-incubation of the enzyme with the inhibitor followed by 
measurement of the initial velocities of substrate hydrolysis as 
a function of pre-incubation time. In this way, increasing the 
concentration of diplacone led to a decrease in both the initial 
velocity (Vi) and the steady-state rate (Vg) (Figs, la, lb and 
Supplementary Fig. S3a). The progress curves obtained using 
various concentrations of the inhibitors were fitted to equation 
(2) (Supplementary Fig. S3a) to determine Vi, Vs and /Cobs- The 
results for diplacone were fitted to the a slow-binding process 
(equation 3; Supplementary Fig. S3a), illustrating a hyperboHc 
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Figure 1 

Inhibition of bacterial sialidase by diplacone. {a) Time-dependent inhibition of Cp-NanI in the presence of diplacone with varying pre-incubation times. 
{b) Pre-incubation time dependence of the fractional velocity of the enzyme-catalyzed reaction in the presence of varying concentrations of diplacone. 
(c) Dependence of /Cobs on the concentration of dipacone. The /Cobs values determined in {b) were fitted to equation (2) in Supplementary Fig. S3(^?). {d) 
Time course of slow-binding inhibition by diplacone. 
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Table 2 

Structures and inhibitory activities (IC50) of flavonoid-based Cp-NanI inhibitors. 



Compoundt 



Structure 



1C50 (m 



Eriodictyol 



3' - O-Methy Idiplacone 



4' - O-Methy Idiplacone 



OCH3 



Mimulone 



Diplacone 




16.4 



14.3 



12.9 



6.1 



0.110 



t lUPAC nomenclature: eriodictyol, (25')-2-(3,4-diliydroxyplienyl)-5,7-diliydroxy-4-cliroman-one; 3'-0-methyldipla- 
cone, 6-[(2£')-3,7-dimetliylocta-2,6-dienyl]-5,7-diliydroxy-2-(3-liydroxy-4-metlioxyplienyl)-2,3-diliydrocliromen-4-one; 
4'-0-methyldiplacone, 6-[(2£')-3,7-dimetliylocta-2,6-dienyl]-5,7-diliydroxy-2-(3-hydroxy-4-methoxyplienyl)-2,3-dihydro- 
chromen-4-one; mimulone, (25')-6-[(2£')-3,7-dimethylocta-2,6-dienyl]-5,7-dihydroxy-2-(4-hydroxyphenyl)-2,3-dihydro- 
chromen-4-one; diplacone, 2-(3,4-dihydroxyphenyl)-6-[(2£)-3,7-dimethylocta-2,6-dienyl]-5,7-dihydroxy-2,3-dihydro- 
chromen-4-one. 



dependence of k^hs on the concentration of bacterial sialidase 
(Fig. Ic). Based on the above kinetic observations, we 
conclude that the reaction follows a time-dependent kinetic 
mechanism (Fig. Id). From the results of the fit, values of 0.201 
± 0.003 and 0.0017 zb 0.004 min"^ were obtained for ks and ke, 
respectively. Using these values derived from kohs gives a 
calculated T^f^P value of 0.065 ib 0.008 which is in 
agreement with the (0.04 jiM) value determined using 
double-reciprocal plots (Lineweaver-Burk and Dixon plots; 
Supplementary Fig. S4). The hyperbolic dependence observed 
with the inhibitors suggests that diplacone inhibits sialidase by 
the rapid formation of an enzyme complex (E-I) which slowly 
isomerizes to form a modified enzyme complex (E*-I) (Fig. 1 
and Supplementary Fig. S3b). 



3.2. Structure of Cp-NanI in complex with diplacone 

Comparison of the structures of sialidases from bacteria 
(Crennell et al, 1993, 1994; Gaskell et al, 1995), viruses 



(Varghese et al., 1983; Burmeister et al., 
1993; CrenneU et al., 2000) and humans 
(Chavas et al., 2005) suggested that 
sialidases share a highly conserved 
active site located within the six-bladed 
y^-propeller fold (Supplementary Table 
S3 and Supplementary Fig. S5). To 
investigate the mode of binding of 
diplacone to the enzyme, we determined 
the structure of Cp-Nanl in complex 
with diplacone. The structure of Cp- 
NanlcD in complex with diplacone 
shows little deviation from previously 
pubHshed structures. The y^-propeller 
fold within the apo structure (PDB 
entry 2vk5; Newstead et al, 2008) and 
the structure of Cp-NanI in complex 
with the substrate Neu5Ac (A^-acetyl- 
neuraminic acid; PDB entry 2bf6; 
Newstead et al. , 2008) show root-mean- 
square deviations (r.m.s.d.s) relative 
to our structure of 0.12 and 0.11 A, 
respectively (Figs. 2a and 3a). Consis- 
tent with our kinetic analysis, diplacone 
was bound to the enzyme active site. 
The principal interactions between the 
inhibitor and enzyme were hydrogen 
bonds and hydrophobic interactions. 
The 3^4^5,7'-hydroxyl-4-oxo group and 
the O atom at the 1 -position of dipla- 
cone form direct or indirect hydrogen 
bonds with side chains of the tri-arginyl 
cluster (Arg266, Arg555 and Arg615) 
as well as residues Arg285, Asp328, 
Phe353, Tyr485, Gln493, Thr538, Tyr655 
and the main chain of Ala292, in addi- 
tion to three water molecules (Wl, W2 
and W3) (Fig. 2h and Supplementary 
Fig. S6). The flavanone backbone of diplacone also displays 
hydrophobic contacts with Ile327, Phe347, Phe460, Tyr485 and 
Tyr655, which orientate the inhibitor correctly and stabilize 
the bound conformation. 

The crystal structure also provides an explanation for the 
weak inhibition shown by the parent eriodictyol and three 
other geranylated flavonoids (3^-(9-methyldiplacone, A'-O- 
methyldiplacone and mimulone) studied in this manuscript 
(Table 2; Fig. 2b). The geranyl group made a suitable hydro- 
phobic interaction with four aromatic residues: Phe353, 
Trp354, Phe460 and Tyr587. Therefore, diplacone (IC50 of 
0.11 \iM) showed an approximately 150-fold higher inhibitory 
effect compared with eriodictyol lacking the geranyl group 
(IC50 of 16.4 \iM). The hydroxy 1 group at C3^ of diplacone also 
made additional hydrogen bonds with the main-chain N atom 
of Ala292. Introducing methylation of this position (3^-(9- 
methyldiplacone, IC50 of 14.3 |iM) or the adjacent C4' (A'-O- 
methyldiplacone, IC50 of 12.9 jiM) would disturb this 
hydrogen bond and cause steric hindrance with the main chain 
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of Ala292 and the side chains of Ile327 and Phe347. On the 
other hand, ehminating the C3' hydroxyl group (mimulone, 
IC50 of 6.1 |iM) does not lead to structural clashes with the 
protein but prevents hydrogen-bond formation. Therefore, 
maintaining the geranyl and Cy hydroxyl group is important 
for the further development of diplacone-based Q?-Nanl 
inhibitors. 

3.3. Comparison with apo and substrate complex of Cp-NanI 

Although the overall structure of Cp-NanlcD in complex 
with diplacone is similar to those of the apo enzyme and the 
complex with NeuSAc, there are interesting structural differ- 
ences between the apo enzyme and the complexes with 
NeuSAc and diplacone that provide useful information 
regarding the development of natural geranylated flavonoid- 



based synthetic inhibitors (Fig. 3b). Asp291 in the diplacone 
and NeuSAc complexes, which is a mechanistically important 
acid/base catalyst (Crennell et al, 1993), rotates by approxi- 
mately 120° compared with the apo structure to make a 
hydrogen bond to the 3^-hydroxyl of diplacone or the 
4^-hydroxyl of NeuSAc. Moreover, the rotation of Asp291 also 
induces a conformational change of Ile327 in the apo enzyme 
compared with both complexes. Relative to the diplacone- 
bound complex, in the NeuS Ac-bound complex Ile327 moves 
to make a suitable hydrophobic interaction with the y^-carbon 
of Asp291 and the S^-acetamide moiety of the NeuSAc 
molecule. Apo Q7-NanIcD and the diplacone complex showed 
a hydrogen-bond interaction mediated by one water molecule 
(Wl), the position of which nearly overlaps with the 
4^-hydroxyl group of NeuSAc. On the other hand, the 
hydrogen bond between Asp328 and the 4^-hydroxyl group of 




Figure 2 

The Q?-Nanl catalytic site with diplacone. {a) Overall structure of the Cp- 
Nanl catalytic domain bound to diplacone. {b) Details of the mode of 
binding. Diplacone and three water molecules (Wl, W2 and W3) are 
shown as green sticks and red spheres, respectively. Hydrogen bonds are 
displayed as dashed lines. 




{a) 



Diplacone 
NeuSAc 




Figure 3 

Structural comparison of apo NanlcD? NanlcD-NeuSAc and NanlcD- 
diplacone. {a) Overall superposed structures of apo NanlcD (orange) and 
the NanlcD-NeuSAc (cyan) and NanlcD-diplacone (yellow) complexes. 
{b) Catalytic site. Diplacone and NeuSAc molecules are displayed as 
green and cyan sticks, respectively. Water molecules (Wl, W2 and W3) 
are shown as red spheres. 
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NeuSAc seems to be responsible for inducing the conforma- 
tional change of Ile327. Q7-NanIcD in complex with diplacone 
has the same side-chain orientation of Asp291 as in the 
substrate (NeuSAc) complex. In addition, an extra water 
molecule (Wl) rather than the 4^-hydroxyl group of NeuSAc 
stabilizes the diplacone complex. Therefore, together with 
the kinetic inhibition assay, the structure of the Q?-NanIcD- 
diplacone complex suggests that diplacone binds to the same 
substrate-binding site and acts as a competitive inhibitor. 

4. Discussion 

We have extracted natural geranylated flavonoids from 
P. tomentosa and shown that they exhibit notable kinetic 
properties against siahdase from C. perfringens, a bacterium 
causing various gastrointestinal diseases. We also found that 
diplacone is the most effective inhibitor among these flavo- 
noids. This inhibitor also shows a time-dependent competitive 
inhibition pattern. The crystal structure of the Cp-NanI cata- 
lytic domain in complex with diplacone rationalized how the 
geranyl and C3' hydroxyl groups of diplacone contribute to 
the stability of the enzyme-inhibitor complex (Table 2). 

From the structural information described above, we 
suggest further modification strategies of diplacone to 
improve its inhibitory properties against C. perfringens Nanl. 
In the Q7-NanI-diplacone complex, two water molecules (Wl 
and W2) mediate hydrogen bonds between Q7-Nanl residues 
(Arg285 and Asp328 with Wl and Gln493, Thr538 and Tyr587 
with W2) of the catalytic site and diplacone atoms (Ol with 
Wl and C4 with W2) (Table 2 and Fig. 2b). Introducing polar 
moieties such as carboxylic and amino groups at the Ol or C4 
positions of diplacone may possibly improve binding affinity in 
the Q7-Nanl catalytic site through the formation of additional 
hydrogen bonds. Addition of hydrophobic groups, including 
phenyl, benzyl and cyclohexyl groups, to C3 of diplacone may 
also increase the binding affinity via hydrophobic contact with 
the hydrophobic patch (Phe353, Trp354 and Phe460) of the 
enzyme. Furthermore, aryl groups (phenyl or benzyl) may 
form a potential dipole-dipole interaction with Phe353. 

Because siahdases are expressed not only in C. perfringens 
but also in human lysosomes (Neul and Neu4), cytosol (Neu2) 
and membranes (Neu3), diplacone could cause undesirable 
effects in humans because it may inhibit human sialidases. To 
predict how diplacone would interact with human silalidase 
counterparts, we modelled human sialidases (Neul, Neu3 and 
Neu4) using the previously determined Neu2 structure and 
compared the specificity of diplacone with C. perfringens and 
human sialidases (Fig. 4a and Supplementary Fig. S2). The 
Q7-Nanl catalytic domain with bound diplacone was super- 
imposed with the crystal structure of human Neu2 with bound 
2-deoxy-2,3-dehydro-A^-acetylneuraminic acid (NeuSAc) and 
models of human Neul, Neu3 and Neu4 that were constructed 
using Neu2 as a template. The geranyl group of diplacone is 
essential for inhibition of the bacterial enzyme as it occupies a 
hydrophobic pocket containing Phe353, Thr487 and Tyr587. In 
human Neu proteins the corresponding residues are polar and 
would not interact favourably with the geranyl group (Fig. 4b). 



Consequently, the structural comparison between the human 
siahdases Neul-Neu4 and the Q7-NanI-diplacone complex 
suggests that the interaction between human siahdases and 
diplacone is likely to be unfavourable owing to polar or ionic 
repulsion. 

On the other hand, despite extensive efforts in drug 
development targeting influenza neuraminidase, a series of 
drug-resistant mutants from influenza viruses have been 
reported, such as the I223R/H275Y dual mutation of neur- 
aminidase from influenza A HlNl, which shows enormous 
resistance towards oseltamivir (approximately 7500 times that 
of the wild type; van der Vries et al, 2012). To overcome drug 
resistance, various natural flavonoid compounds containing a 
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Figure 4 

Structural comparison of Cp-NanlcD, the human Neu2-Neu5Ac2en 
complex structure and homology models of human Neul, Neu3 and 
Neu4. {a) Overall superposed structures of Cp-NanlcD (green), the 
human Neu2-Neu5Ac2en complex structure (grey) and homology 
models of human Neul (purple), Neu3 (brown) and Neu4 (magenta). 
{b) Detailed view. Diplacone and Neu5Ac2en molecules are displayed as 
green and grey sticks, respectively. Water molecules are shown as red 
spheres. 
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flavanone backbone with inhibitory effects on viral neur- 
aminidases (Grienke et al, 2012) may be useful. Because 
diplacone has a flavanone backbone, structural information on 
the Q?-NanI-diplacone complex would help in the discovery 
of viral neuraminidase inhibitors. In this light, we compared 
the structure of the Cp-NanlcD-diplacone complex with two 
structures of drug-resistant mutants of viral neuraminidases 
(the I223R mutant of HlNl and the H274Y mutant of H5N1) 
to examine the use of diplacone against influenza viral neur- 
aminidases. Despite the low sequence identity and low 
homology of viral neuraminidases from influenza A (lA-NA) 
and B (IB-NA) with Q7-NanIcD (identity: lA-NA, 11.56%; 
IB-NA, 12.55%; homology: lA-NA, 29.15%; IB-NA, 29.53%), 
their overall structures are weU conserved (Supplementary 
Figs. SI and S7a). Moreover, diplacone seems to be accessible 
regardless of the I223R or the H275Y mutation. The hydroxyl 
group at the C4 position of diplacone is able to interact with 
Glu276 of the H275Y mutant as in the wild type, which 
suggests that diplacone may be a potential lead compound for 
the development of inhibitors of drug-resistant mutants of 
viral neuraminidases (Supplementary Fig. Sib). 

In summary, we isolated the geranylated flavonoids dipla- 
cone, mimulone and 4^-O-methyldiplacone from P. tomentosa 
and showed their inhibitory activities against Q?-Nanl using 
comparative activity screening. The first crystal structure of 
sialidase with a geranylated flavonoid provides structural 
insights into the binding mode of natural flavonoid-based 
inhibitors at atomic resolution. The higher inhibitory potency 
of diplacone relative to the other compounds studied here can 
be ascribed to two principal factors. Firstly, the 3^-hydroxyl 
group stabilizes the enzyme-inhibitor complex through 
hydrogen bonding to the main chain of Ala292. Secondly, 
additional hydrophobic interactions are observed between the 
geranyl group and the hydrophobic pocket formed by Phe353, 
Trp354 and Tyr587 in the active site. Importantly, in the four 
variants of human sialidase this pocket has a significantly 
different polarity and would not be expected to interact 
favorably with the geranyl group. Glulll, Tyrl79 and TyrlSl 
of human Neu2 are not conserved in bacterial and viral 
siaHdases and are considered as important sites for the design 
of selective drugs targeting bacteria and viruses (Burmeister 
et al, 1993). However, our biochemical and structural data 
suggest that the hydrophobic binding pocket of Cp-NanI 
which binds the geranyl group within diplacone may also be 
exploited as a means to derive selectivity for bacterial versus 
human sialidases. Our results will provide valuable informa- 
tion for the design of new selective antibacterial or antiviral 
agents using flavonoids. 
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